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ABSTRACT 
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The central source of the supernova remnant PKS 1209-52 was observed with 
the Advanced CCD Imaging Spectrometer aboard Chandra X-ray observatory on 
C$ ■ 2000 January 6-7. The use of the Continuos Clocking mode allowed us to perform 

the timing analysis of the data with time resolution of 2.85 ms and to find a period 
P = 0.42412927 ± 2.3 x 10~ 7 s. The detection of this short period proves that the 
CN [ source is a neutron star. It may be either an active pulsar with unfavorably directed 

radio beam or a truly radio-silent neutron star whose X-ray pulsations are caused by a 
J> ' nonuniform distribution of surface temperature. To infer the actual properties of this 

^ ■ neutron star, the period derivative should be measured. 
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1. Introduction 



Radio-quiet compact central objects (CCOs) of supernova remnants (SNRs) have emerged 
recently as a separate class of X-ray sources (see, e.g., Kaspi 2000, for a review). The nature 
^ ■ of these sources, which are expected to be either neutron stars (NSs) or black holes formed in 

supernova explosions, remains enigmatic. They are characterized by soft, apparently thermal, 
X-ray spectra and a lack of visible pulsar activity and optical counterparts. One of the most 
important properties, which potentially allows one to elucidate the nature of these sources, is 
periodicity of their radiation. Some of these sources show periods in the range of 6-12 s, much 
longer than those of typical radio pulsars. These object form a subclass of anomalous X-ray 
pulsars (e.g., Gotthelf & Vasisht 2000), which have been interpreted as magnetars (NSs with 
superstrong magnetic fields — see Thompson 2000). A putative period of 75 ms was proposed 
by Pavlov, Zavlin, & Trimmer (1999) for the Puppis A CCO; however, a low significance of that 
detection requires it to be confirmed by other observations. An unusually long period of 6 hours 
has been reported for the RCW 103 CCO (Garmire et al. 2000a), the origin of this periodicity is 



1 Max-Planck-Institut fiir Extraterrestrische Physik, D-85740 Garching, Germany; zavlin@xray.mpe.mpg.de 

2 The Pennsylvania State University, 525 Davey Lab, University Park, PA 16802, USA; pavlov@astro.psu.edu 



- 2 - 



unclear. No periodicity has been found for the recently discovered CCO of the Cas A SNR (Pavlov 
ct al. 2000; Chakrabarty et al. 2000). 

In this paper we report detection of a period of IE 1207.4-5209, the CCO of PKS 1209-52, a 
barrel-shaped radio, X-ray, and optical SNR (also known as G296. 5+10.0). From the analysis of 
radio and optical observations of this SNR, Roger et al. (1988) estimated its age ~ 7000 yr, with 
an uncertainty of a factor of 3. A recent estimate of the distance to PKS 1209-52, d = 2.1+JJ kpc, 
is given by Giacani et al. (2000). Estimates of the interstellar hydrogen column density from the 
radio, optical, and UV data yield nn,2i = ^_ff/(10 21 cm" 2 ) ~ 1.0-1.8 (see Kellet et al. 1987; Roger 
et al. 1988; Giacani et al. 2000), consistent with a distance d.2 = d/{2 kpc) ~ 1. 

The point source IE 1207.4-5209 was discovered with the Einstein observatory (Helfand & 
Becker 1984), 6' off-center the 81' diameter SNR. Mereghetti, Bignami & Caraveo (1996) and 
Vasisht et al. (1997) showed that the ROSAT and ASCA spectra of IE 1207.4-5209 can be 
interpreted as blackbody emission of T ~ 3 MK from an area with radius R ~ 1.5 c?2 km and 
suggested that this radiation comes from hot spots on NS surface. The spots may be heated either 
by dissipative heating in the NS interior or by the bombardment of polar caps by relativistic 
particles from the NS magnetosphere if IE 1207.4-5209 is an active pulsar. However, the former 
hypothesis can hardly explain the small sizes of the hot spots, even with allowance for large 
anisotropy of thermal conductivity of the magnetized NS crust. The latter heating mechanism is 
also in doubt because of the absence of radio and 7-ray emission from IE 1207.4-5209. 

From observations at 4.8 GHz, Mereghetti et al. (1996) found an upper limit of ~ 0.1 mJy 
on the radio flux from IE 1207.4-5209. They also set a deep limit of V > 25 for an optical 
counterpart in the Einstein HRI error circle, which supports the hypothesis that IE 1207.4-5209 
is indeed an isolated NS. 

Zavlin, Pavlov & Trimmer (1998) reanalyzed the ROSAT and ASCA data, fitting the observed 
spectra with NS atmosphere models. They have shown that the hydrogen atmosphere fits yield 
more realistic parameters of the NS and the intervening hydrogen column than the traditional 
blackbody fit. In particular, for a NS of mass 1.4 M Q and radius 10 km, they obtained a NS 
surface temperature T e g = (1.4-1.9) MK and distance d = 1.6-3.3 kpc, versus T = (4.2-4.6) MK 
and implausibly large d = 11-13 kpc for the blackbody fit, at a 90% confidence level. The 
hydrogen column density inferred from the atmosphere fits, nn,2i = 0.7-2.2, agrees fairly well 
with independent estimates obtained from UV observations of nearby stars, radio data, and X-ray 
spectrum of the shell of the supernova remnant, whereas the blackbody and power-law fits give 
considerably lower and greater values, riH,2i = 0.2-0.4 and nn,2i = 5.2-7.0, respectively. The 
NS surface temperature inferred from the atmosphere fits is consistent with standard NS cooling 
models. 

All the previous observations failed to detect pulsations of X-ray radiation from IE 1207.4- 
5209. An upper limit of 18% (at a 95% confidence level) on flux modulation was set by Mereghetti 
et al. (1996) from the analysis of the ROSAT PSPC data. Due to its higher sensitivity and the 



-3- 



possibility of continuous observations, the Chandra X-ray Observatory is much more capable to 
search for periodicity of X-ray sources. We have employed this capability to search for the period 
of IE 1207.4-5209. The observation and data reduction are described in §2, the timing analysis of 
the data is presented in §3, and some implementations are briefly discussed in §4. 

2. Observation and data reduction 

IE 1207.4-5209 was observed with Chandra on 2000 January 6-7 with the spectroscopic 
array of the Advanced CCD Imaging Spectrometer (ACIS — see Garmire et al. 2000b) in the 
Continuous Clocking (CC) mode. This mode provides the highest time resolution of 2.85 ms 
available with ACIS by means of sacrificing spatial resolution in one dimension. The source was 
imaged on the back-illuminated chip S3. The total duration of the observation was 32.6 ks. Time 
history of detected events reveals that there were three relatively short time intervals with very 
strongly, by an order of magnitude, increased background (background "flares" ) distributed over 
the whole detector. To mitigate the contamination of the source by the flares, we excluded these 
intervals from further analysis, which resulted in the effective exposure of 29,283 s. 

The event arrival times were not properly corrected for the satellite wobbling (dither) by the 
standard pipe-line processing. From our analysis (Sanwal et al. 2000) of Chandra data on PSR 
1055-52, observed in the same mode, we have been aware that the lack of this correction results in 
false periodicities (side peaks in the power spectrum) of a pulsar. Therefore, we corrected the event 
arrival times t in the event file making use of the formula (Glenn Allen, private communication): 

i corr = t + ci [(a- a m ) sinC cos 5 m - (6 - <5 m ) cos (] , (1) 

where c\ = C2 8t, C2 = 3600"/0 / .'4919 = 7318.5605 is the scaling coefficient transforming degrees to 
the detector pixels, 5t = 2.85 ms is the integration time in the CC mode, ( is the roll angle, and 
a — a m and 5 — 5 m (in degrees) are the deviations of right ascension and declination from their 
median values (calculated with the use of the aspect solution file). We have checked that this 
correction removes the artificial periodicities in the data on PSR 1055-52. 

The ID image of IE 1207.4-5209 in "sky pixels" (Fig. 1) is consistent with the ACIS PSF. 
To obtain the source count rate, we extracted 23,337 source+background counts from a ID 
segment of 8 pixel length centered at the source position (10 pixels contain 23,602 counts). The 
background was taken from similar segments adjacent to the ID source aperture. Subtracting the 
background, we find the source countrate of 0.76 ± 0.01 s -1 . A preliminary analysis of the source 
spectrum allows us to estimate the observed source energy flux f x = 2.2 x 10~ 12 erg cm~ 2 s _1 , in 
the 0.1-5.0 keV range. The corresponding (unabsorbed) luminosity is L x « 1.3 x lO 33 ^ erg s _1 , 
and the bolometric luminosity (corrected for the gravitational redshift at the surface of NS with 
M = 1.4M , R = 10 km) is L x « 5 x lO 33 ^ erg s _1 , in agreement with the estimates obtained by 
Zavlin et al. (1998). 



-4 - 



3. Timing analysis 

Previous ROSAT and ASCA observations of IE 1207.4-5209 have not revealed pulsations of 
X-ray radiation from this object. This may be explained by small amount of counts collected in 
those observations. Moreover, those data were spread over very long time intervals, that always 
heavily complicates searching for weak pulsations. 

One of the main advantages of the Chandra observation is that a large number of the source 
counts (about 10 times the number of counts previously detected) has been acquired in a short time 
span. The spectral analysis shows that the background exceeds the source radiation at energies 
above 3 keV. Therefore, for the timing analysis we chose 22,535 counts with energies below 3 keV 
from a 6-pixel segment centered at the IE 1207.4-5209 position. Of these counts, about 98% were 
estimated to belong to the source. To search for pulsations, we used the well-known Z\ (Rayleigh) 
test (Buccheri et al. 1983), which is expected to be optimal to search for smooth pulsations. We 
ran the test in the 0.01-100 Hz frequency range with a step 5f = 3 pHz. The oversampling by 
a factor of 10 (compared to the expected widths of ~ 1/T of the Z\ peaks, where T = 32.6 ks 
is the observational span) was chosen to resolve separate Z\ peaks; it guarantees that the peak 
corresponding to a periodic signal would not be missed. The number of statistically independent 
trials in the chosen frequency range can be estimated as N = 100 Hz x T = 3.26 x 10 6 . The 
test yields only one high peak of Zl mSLX = 65.0 at the frequency = 2.357769 Hz (see Fig. 2). 
Since the variable Z\ has a probability density function equal to that of a x 2 with two degrees of 
freedom, the probability to obtain a noise peak of a given height in one trial is exp(— Zl/2). This 
means that the probability to obtain by chance a peak of Zf = 65.0 in M independent trials is 
p = A/exp(— Z\ max /2) = 2.5 x 10~ 8 , that corresponds to a detection of pulsations at a confidence 
level of C = (1 - p) x 100% = 99.9999975%, or 5.5cr. The next heighest peak is Z\ = 32.0. The 
probability to obtain such a peak by chance in TV trials is 35%. 

To be sure that the pulsations found at the frequency are not associated with instrumental 
effects, we ran the Z\ test at / = /* on the same amounts of counts extracted from a few 
background regions and obtained maximum Z\ equal to 5. We also varied the size of the segment 
for source counts extraction and found that the initial choice of 6-pixel length was optimal to 
produce the maximum Z\ value. We found that contributions from higher harmonics to the value 
of Z 2 are insignificant: Z| = 65.1 and Z| = 65.2. 

To evaluate the pulsation frequency more precisely and find its uncertainty, we employed the 
method suggested by Gregory & Loredo (1996; GL hereafter), based on the Bayesian formalism. 
The method uses the phase-averaged epoch-folding algorithm to calculate a frequency-dependent 
odds ratio, O m (f), which specifies how the data favor a periodic model of a given frequency / 
with m phase bins over the unpulsed model (see Gregory & Loredo 1992 for details of odds ratio 
computations). To weaken the dependence on number of bins, GL suggested to use the odds ratio 
Op C r(f) = X)m=2 X O m (f), with a characteristic number of m max = 10-15. The distribution of 
probability for a signal to be periodic, with frequency / in a chosen frequency range (/i,/2), can 
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be written as 

P(f) 



Oper(/) 



L (i + o; cr )in^ 

where 



(2) 



0; er =(ln|) J^O pcr (f), (3) 

and the probability for the signal to be periodic in (/i, ^2) is 0* cr /(l + 0* cr ). Using PSR 0540-69 
as an example, GL have demonstrated that this method allows one to determine the pulsation 
frequency with much higher accuracy than it can be done with more traditional methods. The 
most probable frequency /o is given by the maximum of O pcr (f) in f2). The uncertainty Sf at 
a given confidence level C = p x 100% can be calculated from the following equation 

rfo+Sf/2 

/ p(f)df = p. (4) 

Jfo-Sf/2 

We implemented the GL method, taking m max = 12 and f 2 -/* = /*- /1 = 200 juHz. The 
frequency dependence of odds ratio is shown in Figure 2. The maximum value, = 6.2 x 10 s , 

is at / = f = 2.3577717 Hz (/ - /* ~ 3 pRz ~ 0.1/T). The uncertainty of f at 68%, 90%, 
and 95% confidence levels is 5f = 1.3, 2.5, and 3.5 /iHz, respectively. This estimate is practically 
independent of choice of frequency range if T f 2 — fi <C /*, where T ~ 3 pRz is a characteristic 
width of the peak of the odds ratio (see Fig. 2). Thus, we finally derive the frequency and the 
period of the detected pulsations 

fo = 2.3577717 ± 1.3 x 10~ 6 Hz, (5) 

P = 0.42412927 ± 2.3 x 10" 7 s , (6) 

at the epoch of 51549.630051303 MJD (TDB). 

The light curve extracted at / = /o (Fig. 3) reveals one broad pulse per period with intrinsic 
source pulsed fraction of f p = 9 ± 2%. Assuming that the detected signal is sinusoidal, we obtain 
an estimate on the pulsed fraction f p = (2Zf mSLX /N) 1 ^ 2 = 7.6% (where TV is the number of counts), 
in a fair agreement with the value calculated from the extracted light curve. Figure 3 indicates 
that the shape of the light curve may vary slightly with photon energy. 



4. Discussion 

The detection of the period, P ~ 424 ms, proves that IE 1207.4-5209, the central compact 
object of the PKS 1209-52 SNR, is the neutron star. 

Let us assume that IE 1207.4-5209 is an active pulsar of an age r within a range of 2-20 
kyr (estimated limits of the SNR age). Then, for a braking index of 2.5, we should expect 
/ ~ //(1.5r) ~ (2.5-25) x 10" 12 Hz s" 1 , P = (4.5-45) x 10~ 13 s s _1 . A crude estimate for the 
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pulsar magnetic field (at the magnetic pole) would be B ~ 6.4 x 10 19 (PP) 1 / 2 ~ (3-9) x 10 13 G, 
close to the critical value B CT = 4.4 x 10 13 G, above which nonlinear QED effects can affect 
the processes in the pulsar magnetosphere. The pulsar spin-down luminosity (rotation energy 
loss) would be E = (2.3-23) x 10 35 /45 erg s _1 . This E may be high enough to power a 
compact synchrotron nebula — for instance, Gaensler et al. (1998) found a radio nebula with a 
radius of 0.3(d/7 kpc) pc around a 100-kyr old PSR B0906-49, whose spin-down luminosity is 
4.9 x 10 35 erg s _1 . A nebula of similar size would have an angular radius of ~ 30" at a distance of 
2 kpc. The ID image of IE 1207.4-5209 puts an upper limit of 3"-4" on the radius of a compact 
nebula around the pulsar. 

If IE 1207.4-5209 is an ordinary radio pulsar whose radio-quiet nature is due to an unfavorable 
orientation of the pulsar beam, the pulsar could be detected in deep radio observations at low 
frequencies, where radio beams are broader. On the other hand, the lack of manifestations of 
pulsar activity may indicate that IE 1207.4-5209 is not an active pulsar, e.g., because a very 
high magnetic field can inhibit the cascade processes in the pulsar's acceleration zone (Baring & 
Harding 1997). In this case, the most natural explanation of the observed pulsations would be 
anisotropy of temperature distribution caused by anisotropic heat conduction in a superstrong 
magnetic field (Greenstein & Hartke 1983). This hypothesis can be verified by a phase-dependent 
spectral analysis which will be done on the same data when the ACIS response is known with 
better precision. Particularly useful for elucidating the nature of the NS would be measuring of 
the period derivative, which requires at least one more observation. For instance, a similar ACIS 
observation in CC mode, taken a year later, would allow us to determine the frequency with an 
accuracy of ~ 2 x 10 -6 Hz and detect a frequency derivative if it exceeds ~ 2 x 10 -13 Hz s _1 , 
which is substantially smaller than the above estimates. 

Our thanks are due to Glenn Allen, who provided the algorithm to correct the event times 
for dither. We are grateful to Gordon Garmire, John Nousek, Leisa Townsley and George Chartas 
for the useful advice on the analysis of ACIS data. This work was partly supported by SAO grant 
GO0-1012X. 
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Fig. 1. — Distribution of detected counts over pixels in ID image in the vicinity of IE 1207.4-5209. 
One sky pixel is equal to / /492 

Fig. 2. — Power spectrum (upper panel) and frequency dependence of odds ratio (lower panel) 
around /* = 2.357769. 

Fig. 3. — Light curves extracted in different energy ranges period Pq. 
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